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Abstract 

We examine the total energy E tot (r) = 2m po i e + Vqcd(^) of the bb system within per- 
turbative QCD up to O(a'g). We extend a previous analysis by incorporating effects of 
the non-zero charm-quark mass in loops. We find that, once the renormalon cancella- 
tion is performed, £tot( r ) agrees well with typical phenomenological potentials for heavy 
quarkonia at distances 0.5 GeV^ 1 < r < 3 GeV -1 . We also examine the perturbative 
predictions for Et t( r ) of other heavy quarkonium systems. Whenever stable predictions 
are obtained, they agree with each other up to an r-independent constant. 



1 Introduction 



Heavy quarkonium systems, such as bottomonium and charmonium, provide an important 
testing ground for theoretical studies on the dynamics of QCD boundstates. For decades, 
theoretical methods for investigating these systems have followed paths different from those 
used for studying the QED boundstates such as positronium or the hydrogen atom. The reason 
is that the non-relativistic boundstate theory based on perturbative QCD, by itself, was not 
successful for describing the heavy quarkonium states. The most serious problem has been the 
fact that perturbative QCD does not reproduce the shape of the static QCD potential VqcdO"), 
which is considered to dictate a dominant part of the dynamics of these system, in the relevant 
region 0.5 GeV -1 < r < 5 GeV" 1 . 

The theoretical method which has been widely used for analyses of the heavy quarkonium 
states is the phenomenological potential-model approach. In this method, one assumes some 
simple form for an effective (non-relativistic) Hamiltonian of the quarkonium system. A phe- 
nomenological potential, which is close to the above QCD potential conceptually, is introduced 
in the Hamiltonian. One parameterises the potential and determines the parameters such that 
various physical observables of the quarkonia are reproduced in this model. It turned out that 
phenomenological potentials determined in this way have more or less similar shapes in the re- 
gion 0.5 GeV" 1 < r < 5 GeV" 1 , which may be represented typically by a Coulomb-plus-linear 
potential. See e.g. Ref. |IJ for a recent analysis based on potential models. 

On the other hand, theoretical approaches, which have foundations on first principles, 
have been developed. Within the framework of Non-Relativistic QCD (NRQCD) or potential- 
NRQCD, the Lagrangian of a heavy quarkonium system is given in a series expansion of some 
small parameter (typically the quark velocity v ~ a s ). The Wilson coefficients of the effective 
theory are determined by matching the theory to the full QCD theory either perturbatively or 
non-perturbatively. The QCD potential, as well as other (sub-leading) potentials which enter 
the effective theory, are determined from the results of lattice calculations or from model cal- 
culations; see p| and references therein. The lattice results, being first principle calculations, 
turn out to be consistent with phenomenologically determined potentials, but at the present 
stage they are still not precise enough. 

Computations of the QCD potential within perturbative QCD also made progress over time. 
The full two-loop corrections with massless quark loops |J, as well as non-zero mass effects in 
quark loops up to the same order ||-@], have been computed. The perturbative expansion at 
r > 0.2 GeV -1 revealed to be very poorly convergent, and also its shape deviates qualitatively 
from an expected Coulomb-plus-linear form in the relevant range. The poor convergence is 
considered to reflect a non-trivial structure of the QCD vacuum. Also, within the context of 
perturbative QCD, this behaviour has been understood using the renormalon language 0. 

Recently, there have been significant developments in the non-relativistic boundstate theory 
based on perturbative QCD. Thanks to the calculations of higher-order corrections @-l L| and 
the understanding of the leading renormalon cancellation [|T^,|T^|, it became possible to predict 
accurately the physical observables of the heavy quarkonium states (in particular, the bottomo- 
nium states) within perturbative QCD |5|,p!0|, |l4HT9j . An essential feature is that the prediction 
for the QCD potential, after incorporating these developments, has become accurate, and that 
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it has reproduced a realistic shape of the potential in the relevant range. Ref. [17] investigated 
this feature in particular. There, the following aspects have been shown: 



(1) When the leading renormalon cancellation is incorporated, convergence property of the 
total energy of the bottomonium system E tot {r) = 2m h vo \ c + VqcdO") improves drastically. 

(2) For simplicity, two hypothetical cases, m c — > and m c — > rrib, have been examined; a 
reliable theoretical prediction for E tot (r) is obtained at r < 3 GeV -1 , and the prediction 
agrees with typical phenomenological potentials in the range 0.5 GeV" 1 < r < 3 GeV" 1 
within theoretical uncertainties. 

(3) The qualitative behaviour of E tot (r) in the above range can be understood as originating 
from an increase of the interquark force due to the running of the coupling constant. 



In this paper we extend the analysis of [17|. First, we incorporate the realistic value of the 



charm mass in the calculation of E tot (r). Since we are interested in the range of r not very 
different from the charm- mass scale 1 / m c , we should properly take into account the dependence 
on m c which enters through loop corrections. In this way we can compare predictions of 
perturbative QCD in the realistic case with phenomenological potentials. Secondly, we compute 
the total energies for quark-antiquark systems other than the bb system and compare them. 
In [JT7J the part of E tot (r) independent of the external quark masses has been examined using 
the interquark force. Here we examine it in a different way. Our purpose is to compare the 
QCD potential (or the corresponding potential in the effective Hamiltonian of the quarkonium 
system), which has been studied in various approaches so far, with the prediction of perturbative 
QCD. We anticipate that, by combining our results with the conventional studies, we would 
be able to obtain a better understanding on the dynamics of the heavy quarkonia from first 
principles. 

In Sec. 2 we compute the total energy of the bb system, incorporating non-zero charm mass 
effects: in Sec. 2.1 we set up our formulas; in Sec. 2.2 numerical analyses are given; in Sec. 2.3 
we discuss uncertainties of our predictions. In Sec. 3 the total energies of other systems are 
examined. Conclusions are given in Sec. 4. Appendices collect some formulas. 



2 Total Energy of the bb System 
2.1 Definitions 

In the MS scheme, it is appropriate to compute the total energy of the bb system in the theory 
which contains 5 flavours. We are, however, interested in the total energy when the distance 
between b and b is much larger than their Compton wavelength, r ^> l/m&. Therefore, we will 
rewrite the total energy in terms of the 4-flavour coupling acg (//) in order to realize decoupling 
of the 6-quark to all orders. 

When we neglect the charm quark mass, the total energy is given by 

Et, mc =o(r) = 2m fe , pole + V QCDA (r)- (1) 
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The relation between the pole mass and the MS mass has been computed up to 3 loops in a 
full theory, which contains heavy flavours and n/ massless flavours [p0|| . (The same relation 
was obtained numerically in |2l| in a certain approximation.) Setting rih = 1 and rewriting the 
relation in terms of the coupling of the theory with n\ = 4 massless flavours only, we findF] 



_ , , 4 af{m b ) 
m btPole = m b < 1 + h 

O 7T 



( 4 ) I— \ 



7T 



(2) 



,MS^MS^ 



;/' (;// ' ) denotes the MS mass renormalised at the MS-mass scale. The QCD 



where rrii 

potential of the theory with n\ massless flavours onlyQ is given, up to (9(a|), by 

4 4 ni V) 



^QCD,n ; (r) 



(ni) i \ . 

I (^-Mj(2^ + a^ ) ) 
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where £ = log(/ir) + 7 B , and $ denote the coefficients of the beta function 
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The constants a\ and dj are given in Appendix A. 

When we include the effects of the non-zero charm-quark mass, the above formula is modified 
as follows: 



Et t {r) 



E 



bb i 

tot,m c =0 \ 



2 5m b)VO \ e + 5Vq C d(^)- 



(5) 



At 0{or s ) the non-zero charm mass correction to the pole mass, 5m b:PO \ e , reads |22 

2 
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+(i - o(i - a Li 2 (o - - io g 2 (o + io g (o iog(i - o - 
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where £ = m c /m b . Its leading term in the limit m c — > (linear approximation) is given by 

(4 4) (m 6 )) 2 



5m 



[2] 

6, pole 



-77V 



(6) 



(7) 



m c ^r0 Q 

* This relation coincides with Eq.(14) of |20| ], which is given numerically (indirectly through (3q and 
Note that, in the other formulas of ^Oj, the coupling of the full theory is used. 

' Due to the decoupling theorem, the perturbative QCD potential of the theory which contains one heavy 
flavour (with mass m) and n; massless flavours coincides with the potential in Eq. (fj) up to C(a|) if we count 
1/r = 0(asm) and if we rewrite the coupling by the coupling of the theory with ni massless flavours only. 
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At C(a|), the complete expression of Sm btPO i e is not known; it has been computed only in the 
linear approximation [H: 



5m 



[3] 

ft, pole 



71 
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where / 2 = 0.470 ± 0.005, b 2 = 1.120 ± 0.010, f x = (log A - log 6 2 )/(log f 2 - \ogb 2 ), h = 
(log A - log/ 2 )/(log& 2 - log/2) and log A = 161/228 + 13C 3 /19 - log 2. In § it has been 



argued that the use of the linear approximation in both Smfj le and &m b ^ po[e is a slightly better 

f2l 

approximation of the full result than to use the exact 5m b le and the linear approximation of 



[3] 

5m [ b l polc . It has been conjectured that the former approximation accounts for the full correction 
121 [31 

1 x - l0i with about 10% accuracy, while the latter approximation accounts for the 



^ m ft,pole + ^ m ft,polc 

full result with about 20% accuracy. In [19|, the theoretical predictions for the bottomonium 
energy levels turned out to be more stable when we used (^^Liel^-^) + (^ m ft 3 poiclm c -»o) as 
compared to 6mfl ole + (5m[ 3 pole | mc _^ ). 

At 0(a 2 s ) the non-zero charm mass correction to the QCD potential, <5Vqcd, is given in 
one-parameter integral form as 



4 a 
3 ~ 



4%) 



V 3tt 



dx f(x) 
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(lot 



(m c r) + 7e + - 
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(9) 



with 
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\fxF 



(10) 



At 0(ag), the correction was computed first in momentum-space in 0. The Fourier transform 
was performed| and SVqq D (t) was obtained also in one-parameter integral form in |5|,|6|]. Both of 
the latter references contain misprints, however; for completeness, we give a corrected formula 
for SVQQ. D (r) in Appendix B]§] 

Since the renormalon cancellation at each order of the perturbative expansion is realized 
only when we use the same coupling constant in expanding m bjPO \ e and VQCD,4( r ), we rewrite 
asijn b ) in terms of as(/J>) using the renormalizat ion-group evolution of the coupling constant: 



(n,)/_ \ (n,)/ \ Ji , Po" n ( V 

a' s '{mi) = ay{fi)<l + — — log — 
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7T 



A 



4 Vmj/ £ 
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* The coordinate-space potential was studied using different integral representations in 
§ The corrected formula has been acknowledged by the author of ||. 
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Furthermore, we also examine the total energy after re-expressing it in terms of the 3-flavour 
coupling and compare it with the 4-flavour coupling case. This is because we are interested in 



the total energy in the range 0.5 GeV 
as well. We insert the relation [f23 



-i < 



^ 5 GeV 1 , where the charm quark may decouple 



q4 4 V) 
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(12) 



into the above E^ t (r) and re-expand in (//). Thus, we will examine the series expansion of 



£; t ft b t (r;m b ,m c ,4 n!) (A«)) in ctp'(fj) up to 0((a^) 3 ) for m = 3 and 4. 

The obtained total energy depends on the scale /i due to truncation of the series at a finite 
order. Following the prescriptions of JI7 |, we will fix the scale fi in the two different ways 
described below: 



(™i)\3^ 



1. We fix the scale fi = /ii(r) by demanding stability of E tot (r) against variation of the scale: 



V-j^E tot (r;rni,as(iJ,)) 



(13) 



2. We fix the scale fi = ^{r) on the minimum of the absolute value of the last known term 
[0(a|) term] of E tot (r): 



d_ 



0. 



(14) 



/i=/i2(r) 



2.2 Numerical Analyses 



In this subsection we take the input value for the coupling constant as ct§ (Mz) = 0.1181 ± 
0.0020 [0]. We evolve the coupling by solving the 3-loop renormalization-group equation 
numerically and match it to the 4- and 3-flavour couplings successively through the matching 
condition [Q.Q For the bottom- and charm-quark masses, we use the values nib = 4.190^9 GeV 
19| and m c = 1.243 GeV []16| |, respectively. (For simplicity we do not change the value of m c as 



a function of otg (Mz)- This is justified, since the dependence of E tot (r) on it is much smaller 
than other theoretical uncertainties; see Sec. |2l 
the linear approximation in this subsection. 

First we examine the scale dependences of E tot and 



We compute both Sm b 2 ^ pole and Sm l b pole m 



.[3] 



ES t (r;mi,a s (ii))\ in Eqs. (|13|) and 



[3] 



(|T^). In the first scale- fixing prescription, the minimal sensitivity scale fii(r) exists only in the 
range r < 3 GeV -1 . In the second prescription, /i = fi2{ r ), the minimum value of \E, 

» The matching scales are taken as rn\, and m c , respectively. 
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Figure 1: Determination of the scale: Total energy and a^-term of the total energy, the scales 
determined with the respective prescription are marked with a short vertical line. (We do not 
take the scales which are located close to the infrared singularity of a s (/i).) The distance r is 
fixed to 2 GeV" 1 ; af(M z ) = 0.1181. 

zero in the range r < 3 GeV -1 , whereas | -E^tot C 7 ") I > in the range r > 3 GeV" 1 . These features 
indicate an instability of the perturbative prediction for E tot (r) at r > 3 GeV™ 1 . Qualitative 
features of the scale dependences are similar for the expansions in the 3-flavour and the 4-flavour 
couplings. (They are also similar to the results of the analysis for m c — > or m c —>■ ||17|| .) 
The range where the prediction is stable extends to slightly longer distances for the 3-flavour 
case, in accord with our naive expectation. The scale-dependence is demonstrated in Fig. [l 
for r = 2 GeV™ 1 and for the expansion in the 3-flavour coupling. We compare the total 
energy determined in the four different prescriptions (scale-fixing prescriptions 1 and 2, and 
expansions in org and org ) in Fig. |2|. All four curves agree well in the range where both scale- 
fixing prescriptions exist. The numerical values of the total energy and the scale are listed in 
Table [1]. We observe a strong cancellation of the leading renormalon between the pole mass 
and the QCD potential: the expansion of E tot (r) is much more convergent than the expansions 
of the individual terms. This is demonstrated in Table 0. 

Since the new aspect of our analysis is the inclusion of the charm-mass effects, we compare 
our result with those of the analysis done in the two hypothetical cases m c — > and m c — > m^. 
This is shown in Fig. [| Our curve lies between those of the two hypothetical cases, which is 
consistent with our expectation. 

Now we compare the total energy E to t{r), as obtained above, with typical phenomenological 
potentials in the literature: 

• A Coulomb-plus-linear potential (Cornell potential) |25| : 




(15) 
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1 2 3 4 5 

r [GeV^ 1 ] 

Figure 2: Total energy for 3-flavour (// = solid, = /i 2 : dashed) and 4-flavour (/jl = fix'- 
dotted, /i = /i 2 : dash-dotted) case. For the 3- (4-) flavour case prescription 1 breaks down at 
r « 3.2(2.75) GeV -1 . Where both scale prescriptions exist, the curves coincide. a${Mz) = 
0.1181. 



r [GeV" 1 ] 


E tot (r) [GeV] 


Mr) [GeV] 


0.5 


8.75 


1.99 


1.0 


9.25 


1.52 


1.5 


9.49 


1.28 


2.0 


9.66 


1.12 


2.5 


9.79 


0.98 


3.0 


9.92 


0.79 



Table 1: Comparison of the total energies and scales for the bb system (3-flavour prescription, 
fi = fi 2 , af{M z ) = 0.1181). 
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n 


M 

m 6,pole 


^ m l,pole 


y[n\ 


V QCD 







4.190 











8.380 


1 


0.657 





-0.493 





0.822 


2 


0.175 


0.028 


-0.298 


-0.032 


0.076 


3 


0.182 


0.044 


-0.391 


-0.084 


-0.024 



Table 2: Renormalon cancellation in the bb system (r = 1 GeV~ , /i = /ii, 3-flavour prescrip- 
tion, Os\m z ) = 0.1181): The a^-expansion of E tot has much better convergence properties 
than the expansions of the individual terms. represents the order term of X. 

^6, P oie and Vqcd do not contain the m c -effects, these are given as 5m fe pole and SVqcd- The 
units are GeV. 




1 2 3 4 5 

r [Ge^ 1 ] 

Figure 3: Charm mass effects: Comparison of the QCD potential for realistic m c with the two 
limiting cases m c — ■> (1 heavy flavour, 4 light flavours) and m c — * mb (2 heavy flavours, 3 
light flavours) for af(M z ) = 0.1181. 
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with k = 0.52 and a = 2.34 GeV- 1 . 
A power-law potential p6| : 

V{r) = -8.064 GeV + (6.898 GeV)(r x 1 GeV) 01 . (16) 

A logarithmic potential |27 |: 



V(r) = -0.6635 GeV + (0.733 GeV) log(r x 1 GeV). (17) 

In Fig. [| we compare these potentials with the perturbative computation of E tot (r) for the 
input values a^(Mz) = 0.1181 ± 0.0020. An arbitrary constant is added to each potential, as 
well as to E tot (r), such that it coincides with E tot (r) for ag(M z ) = 0.1181 at r = 1 GeV _1 |j] 
This is because the constant part (r-independent part) of each phenomenological potential 
is not determined well. As can be seen, E tot (r) and the phenomenological potentials agree 
inside the errors estimated from the next-to-leading renormalons ±^Aqcd(Aqcd ■ f) 2 (taking 
Aqcd = 300 MeV, indicated by error-bars), and the agreement is better for a larger value of 
the input parameter a^\M z )- 

In order to quantify the differences between the phenomenological potentials and the QCD 
potential for different values of as(Mz), we define a weighted difference between E tot (r) and a 
potential V(r) as 

A„,[£ tot (r) - V{r)\ = min f ' drr n \E tot (r) - V{r) - of. (18) 

The minimum value of the integral is taken as we vary an arbitrary constant added to the 
potential. We examined the difference A n between E tot (r) and each of the phenomenological 
potentials, as well as A n between E tot (r) and the average of the three phenomenological po- 
tentials, for n = —1, and +1. We found similar qualitative features: A n becomes small for a 
larger value of (Mz), inside the present error, around 0.1191-0.1201. We demonstrate this 
in Fig. A n=0 with r = 1 GeV -1 and r\ = 4 GeV -1 is shown for each of the potentials and 
for the averaged potential, varying the value of as(M z ). 



2.3 Error Estimates 

There are several uncertainties in our theoretical prediction for the total energy E^ t (r). 

(i) Errors of the input parameters ag{Mz), rh~b andm c : 

In Fig. |] the dependence of E tot {r) on ces{M z ) is shown. As as{M z ) increases, -E"t t(r) 
becomes steeper due to an increase of the interquark force |17|]. Variation of E tot at 

" The choice of r, where the potentials are made to coincide, is to some extent arbitrary. (Our choice 
r = 1 GeV -1 corresponds approximately to the root-mean-square radius of the T(IS') state, the heaviest state 
used to determine the phenomenological potentials.) The important point here is that we can indeed choose 
an additive constant such that Etot(f) can be made to agree with the phenomenological potentials inside the 
estimated errors. 
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1 2 3 4 5 

r [GeV- 1 ] 

Figure 4: Comparison of the QCD potential to various models. The QCD potential is given 
for three values of a^\M z ), constants are added to make all curves coincide with the QCD 
potential (af(M z ) = 0.1181) at r = 1 GeV" 1 . 



0.007 
0.006 
. 005 
0.004 
0.003 
0.002 
0.001 

af{M z 




0.1181 0.1191 0.1201 0.1211 



Figure 5: A comparison of the QCD potential for different ctg(Mz) to various models. The 
height of the bars indicate the minimised difference A n= o between the QCD potential and the 
model when varying the constant term. The four bars in each group correspond to (from left 
to right): Cornell potential, power-law potential, log potential and the average of the former 
three. 



10 



r = 3 GeV -1 is about ±90 MeV with respect to its value at r = 1 GeV" 1 , corresponding 
to a variation cv§\Mz) = 0.1181 ± 0.0020. The present uncertainties of the 6-quark and 
c-quark MS masses are about ±30 MeV |L9| and ±100 MeV [16|,|28|, respectively, and also 



their central values are correlated with the value of as(M z ). The dependence of E tot (r) 
on Mb for small variations of rn~b is practically a shift of E tot (r) by 2 Arab (independent 
of r). We have confirmed this feature up to a larger variation of Am), = ±60 MeV, 
corresponding to more conservative error estimates in the literature. The dependence of 
Etot{ r ) on m c is weak and negligible (^^^(r)! < 10 MeV for Am c = ±100 MeV), since 
the charm mass effects are 0(a 2 s ) and beyond. Hence, of these only the dependence on 
a,s{Mz) matters in the comparison of E'tot(^) and the phenomenological potentials. 

(ii) Unknown higher-order corrections: 

We may estimate the higher-order uncertainties from the size of the next-to-leading renor- 
malon |7|[l7|. We show its typical size ±|Aq CD (Aq CD • r) 2 by error-bars in Fig. f| taking 
Aqcd = 300 MeV. We have also checked these estimates by computing the 0(a A s ) correc- 
tion to E tot (r) in the large-/?o approximation. The corrections from the (9(a|)-term are 
7, 47, 68 MeV at r = 1, 2, 3 GeV -1 , respectively. 



(iii) Approximations in 5mb t 



pole- 



In Sec. ^]2| the charm mass corrections to the pole mass, Sm b 2 ^ pole and Smf^ pole , have 
been computed in the linear approximation. We estimate uncertainties induced by this 
approximation by using the full correction Eq. (||) for 5m^ po[e instead, while keeping 

8m bpole in the linear approximation. We find that E tot (r) varies by —7, —10, —17 MeV 
at r = 1, 2, 3 GeV -1 , respectively. 

(iv) Non-perturbative effects: 

Up to now there exists no reliable way to estimate entire non-perturbative effects on the 
QCD potential accurately 0. Moreover, in principle, the size of non-perturbative effects 
will depend on the specific perturbative scheme used in the computation of the QCD 
potential or the total energy. Here, we do not attempt to compute non-perturbative effects 
separately, but rather consider the differences between our predictions (with respective 
perturbative errors) and the phenomenological potentials as orders of magnitude estimates 
of non-perturbative effects in our computational scheme. From our present results, we may 
consider that non-perturbative effects would be absorbed into perturbative uncertainties 



of our predictions. This is consistent with the observation of [19 



Besides, when we compare our theoretical predictions with phenomenological potentials, we 
should take into account the following uncertainty: 



Phenomenological Potentials vs. E^ t (r): 

Our ultimate goal would be to compare the exact QCD potential with the perturbative 
prediction (after subtracting renormalons) . The phenomenological potentials are not di- 
rect physical observables but are determined under certain model assumptions. Since 
separation of the QCD potential from the rest of the interactions (e.g. 1/r 2 potential, 
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relativistic corrections) is not absolutely clear in the phenomenological approaches, there 
may be substantial corrections between the phenomenological potentials and the QCD 
potential or the total energy E^ t (r). The relation can only be clarified by detailed com- 
parisons of the predictions of perturbative QCD with the experimental data for various 
physical observables of the bottomonium states, such as the energy spectrum, decay rates, 



level-transition rates, etc. (Studies on the energy levels have been initiated in [16. [I 
These detailed comparisons, however, are beyond the scope of the present paper. 



3 Other Systems 

Let us also examine the total energies of other quark-antiquark systems. Theoretically we 
may expect that a limit of the potential energy exists when we send the masses of quark and 
antiquark to infinity. Empirically both the bottomonium and charmonium spectra can be 
reproduced well with the same phenomenological potential. Thus, it would be interesting to 
examine whether the total energies for different systems coincide up to an additive constant. 
For the be system, we consider 

E^ t (r) = m b:PO i c + <5m fe>po i e + m CjPO i e + VQ CD>4 (r) + <5VqcdO), ( 19 ) 

where 



^■c,pole 




(20) 



We examine the expansions of E^ t (r) both in the 4-flavour and 3-flavour couplings, and in the 
two scale-fixing prescriptions Eqs. ( |i~3f ) and (|i~4l) . We find that, at distances 0.5 GeV -1 ^ r 



< 1.8 GeV -1 , /Ui(r) exists in the first prescription and the minimim value of \E^j t \ is zero in 
the second prescription, indicating that stable theoretical predictions for the total energy are 
obtained in this region. A comparison in Fig. |] shows that the predictions for the total energies 
of the be system and bb system agree within the uncertainties of the predictions. We also note 
that, even at r > 1.8 GeV -1 , where only the second scale-fixing prescription exists, the curves 
for the be system in the 4-flavour and 3-flavour couplings agree with each other and they are 
also consistent with that of the bb system. See also Tables ^| and £| for numerical values. 
Furthermore we examined the total energy of the cc system, 

E*(r) = 2m c , po]B + V QpD ,3(r), (21) 

in a similar way. We, however, obtained stable predictions only in the very narrow vicinity of 
r ~ 1 GeV -1 and consequently we could not compare the shape of E tot {r) with that of the 
bb or be system. We believe that this is caused by the typical scales for the cc system being 
too low to give reliable perturbative expansions. (We may compare this feature with that of 
the perturbative calculation of the energy levels of charmonium ||16|| : only the IS* levels can be 
computed reliably.) 
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r [GeV -1 ] 


E tot {r) [GeV] 


// 2 (r) [GeV] 


0.5 


5.35 


4.08 


1.0 


5.85 


1.31 


1.5 


6.12 


0.97 


2.0 


6.34 


0.76 


2.5 


6.48 


0.83 


3.0 


6.55 


0.93 



Table 3: Comparison of the total energies and scales for the be system (3-flavour prescription, 
jjL2). The predictions are stable in the range 0.5 GeV -1 ^ r % 1.8 GeV -1 . 



n 


M 

m 6,polc 




M 

m c,pole 


K QCD 


K QCD 


^tot 





4.190 





1.243 








5.433 


1 


0.730 





0.217 


-0.547 





0.399 


2 


0.153 


0.035 


0.211 


-0.321 


-0.039 


0.038 


3 


0.203 


0.054 


0.304 


-0.466 


-0.109 


-0.014 



Table 4: Renormalon cancellation in the be system (r = 1 GeV , \i = fj>i, 3-flavour prescrip- 
tion). Notations are same as in Table [| 




12 3 4 

r [GeV -1 ] 



Figure 6: Total energy of the S c -system for 3-flavour (/i = n\\ solid, \i = fi2- dashed) and 
4- flavour (// = dotted, // = /i 2 : dash-dotted) prescription. For comparison the energy 
of the T-system (3 flavours) is shown (shifted so that it coincides with the B c energy at 
r = 1 GeV -1 ). For all curves a^\M z ) = 0.1181. Error-bars indicate ±|Aq CD (Aq CD ■ r) 2 with 
A QC D = 300 MeV 



13 



Conversely we examined the total energy of the bb system as we increase the bottom quark 
mass artificially. When we do this, we find that the region of r, where we can make a stable 
prediction of E to t{r), shifts to shorter distances. In other words, we can still make a stable 
prediction of the shape of the potential from the distance r ~ 1 /m& up to the region which is 
relevant to the formation of the (hypothetically heavy) bottomonium states. This is reasonable, 
since a heavier quarkonium has a smaller radius, and the prediction should be stable for this 
smaller radius. On the other hand, we do not always have stable predictions for E tot (r) in 
the region of our interest, 0.5 GeV -1 < r < 5 GeV -1 . We may understand this property as a 
manifestation of a multi-scale problem. The pole mass contains powers of log(/i/mb), while the 
QCD potential contains powers of log(/ir). When and r _1 are very different, it generally 
becomes more difficult to cancel two different types of logarithms, so it becomes difficult to find 
a scale /i which stabilizes the theoretical prediction. For a (hypothetical) value of below 
about 23 GeV, we still have an overlap between the region of stable predictions for the total 
energy and the above region of our interest. Where both scale-fixing prescriptions exist, the 

To] 

predictions coincide. (The distance r where the minimum of | -E"tot I deviates from zero in the 
second prescription is close to the distance beyond which the first scale /ii(r) does not exist.) 
In the region where both prescriptions exist, the predictions agree with E^ t (r) computed with 
the realistic 6-quark mass, up to an additive constant .Q At larger r the larger mass produces a 
smaller total energy. This tendency is expected from the higher-order analysis in the large-/3o 



approximation [17 . Also this agrees qualitatively with our results for the be against bb system: 



from the point where the curves split, E^ t is larger than E^ t , although the curves later cross 
(Fig. P). These features are demonstrated in Fig. [7|, where the total energy for m& = 15 GeV is 
shown and compared with the one with the realistic 6-quark mass. 



4 Conclusions 

We have analysed the total energy of the bb system incorporating the non-zero charm-quark 
mass effects. We observed an improvement of convergence of the perturbative expansion, once 
we perform the cancellation of the leading renormalons; we obtained stable theoretical predic- 
tions for the total energy E^ t (r) at r < 3 GeV -1 . These features are qualitatively the same as 
those observed in |T7]]. We compared the total energy with typical phenomenological potentials. 



They agree in the range 0.5 GeV 1 ^ r < 3 GeV 1 within the estimated theoretical uncertain- 
ties. The agreement becomes very good when the value of the input parameter ag\Mz) is 
large (inside the error bands), 0.1191-0.1201. From these results, we may conclude that, in 
analysing the nature of the bottomonium states, using the perturbative prediction for E^ t (r) 
is, at least, as good as using phenomenological potentials. We may compare our result with 
the comprehensive analysis of the bottomonium spectrum which includes full corrections up to 
0(l/c 2 ) ||19|| . There, a smaller value around 0.1161 for a^(Mz) was favored. Therefore, we 



see that the interactions other than E to t(r) play non-negligible roles for the predictions of the 

* Of course, the agreement is better at shorter distances, where the perturbative predictions tend to be more 
reliable. At r < 1.5 GeV -1 the expansion in the 4-flavour coupling tends to be more stable than the expansion 
in the 3-flavour coupling. 
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12 3 4 

r [GeV -1 ] 



Figure 7: Total energy of the bound state T of a hypothetical heavy 6-quark with a mass of 
15 GeV in 3-flavour (// = /ii [breaks down at r « 0.6 GeV -1 ]: solid, /x = /J2: dashed) and four 
flavour (/i = /ii [breaks down at r « 0.8 GeV -1 ]: dotted, /i = /i2: dash-dotted) prescription. 
Where both scale prescriptions exist, the curves coincide. For comparison the energy for the 
T-system (3 flavours) is shown (shifted so that it coincides with the T energy at r = 0.5 GeV -1 . 
For all curves af{M z ) = 0.1181. 
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bottomonium energy levels, with respect to the present theoretical accuracy of perturbative 
QCD. 

We also examined the perturbative predictions for E tot {r) of other systems. For the be 
system, stable predictions are obtained (at least) in the range r < 1.8 GeV" 1 . In this range, 
the total energy agrees with that of the bb system inside theoretical uncertainties. For the cc 
system, we could barely obtain stable theoretical predictions for the shape of E tot (r), since 
the relevant scales are very low. We also found that, for a heavier (hypothetical) quarkonium 
system, the range where E tot (r) can be predicted reliably, shifts to shorter distances. At any 
event, whenever we may obtain stable theoretical predictions for E to t(r) at r > 0.5 GeV^ 1 , the 
predictions agree with the phenomenological potentials within present theoretical uncertainties. 



Appendix 

A: Parameters 

The constants used in Eqs. (0) and (13) are given by || 
(n) 31 10 

a i - K zz ) 

(n t ) 4343 2 9vr 4 / 1229 52 ( 3 \ 100 2 

a\ u = h 36 7r + 66 Cs ni H nf, 

2 18 S 4 V 27 3 J 81 11 



and 



(23) 



1 32 3 9 6 V 144 18/ 

~ 13.4434- 1.04137 n h (24) 



8462917 652841 vr 2 695 vr 4 575 vr 2 log 2 



4 n;) = + 



93312 38880 7776 162 

22 vr 2 log 2 2 55 log 4 2 220Li 4 (|) : 58 Cs 1439 n 2 ( 3 1975 Cs 



81 162 27 27 432 216 

/ 231847 991tt 2 61tt 4 ll7r 2 log2 27r 2 log 2 2 log 4 2 

+Ul \ 23328 648 + 1944 81 + 81 + 81 

8Li 4 (|) 241C 3 \ 2 / 2353 13tt 2 7( 3 \ 



^ 27 72 ) +ni ^23328 ^ 324 " r 54 ) 

190.391 - 26.6551 m + 0.652691 n]. (25) 
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After correcting misprints in the charm-mass correction to the QCD potential at 0(a 

readsf] 

A (4)/ \ / (4), 



r V 3tt / 



3 r 

{ ~ 2^ dX f<K ^ ( l0g ^J^~ ~ Ei ^ 2mcrx ) _ e 4 " cra Ei(-2m c rx)) - a[ 

c (4) 2 1 

+3 (log(m c r) + 7s + |) + + ^ ^ J 

j^" dz /(x) er^ rx + x /(x) log 



x — vx 



X + \J x 1 — 1 

dx/(x) e - 2mcra (log4x 2 - Ei(2m c rx) - e 4mcra Ei(-2m c rx; 
-(log(m c r) + 7s + - ^ | 

r 57 161 13 

| — (A T(0, 2/ 2 m c r) + h T(0, 26 2 m c r) + log(m c r) + lE + — + — ( 3 



J°° dxf(x) e" 2mcr:r (-8m c rx) +4 



(26) 
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